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Abstract 

Doubly charged Higgs bosons (H^) are a distinctive signature of the Higgs Triplet Model 
of neutrino mass generation. If i? ±± is relatively light (m H ±± < 400 GeV) it will be 
produced copiously at the LHC, which could enable precise measurements of the branching 
ratios of the decay channels i? ±± — > ^ I - . Such branching ratios are determined solely by 
the neutrino mass matrix which allows the model to be tested at the LHC. We quantify the 
dependence of the leptonic branching ratios on the absolute neutrino mass and Majorana 
phases, and present the permitted values for the channels e^e^, e ± // ± and // ± At ± . It is 
shown that precise measurements of these three branching ratios are sufficient to extract 
information on the neutrino mass spectrum and probe the presence of CP violation from 
Majorana phases. 
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1 Introduction 



The firm evidence from a variety of experiments that neutrinos oscillate and possess a small mass 
below the eV scale necessitates physics beyond the Standard Model (SM). Consequently models 
of neutrino mass generation which can be probed at present and forthcoming experiments are 
of great phenomenological interest. In particular, those models which can provide a distinctive 
experimental signature, such as a New Physics particle with a mass of the TeV scale or less, 
are especially appealing in light of the approaching commencement of the CERN Large Hadron 
Collider (LHC). 

Doubly charged Higgs bosons (-ff ±± ) arise in a variety of models of neutrino mass generation 
as members of / = 1, Y — 2 scalar triplets [1-11] and I — 0, Y — 4 scalar singlets [12]. Such 
particles can be relatively light (i.e., with masses of the electroweak scale) and have impressive 
discovery potential at hadron colliders due to their low background signature H ±±L — > if if 
and sizeable cross-sections. The ongoing searches at the Fermilab Tevatron [13, 14] anticipate 
sensitivity to m#±± < 250 GeV for the decay channel if ±± — > if if = e, fx) with the 
expected final integrated luminosities of up to 8 fb _1 . LHC simulations [15, 16] show that 
discovery for m H ±± < 1 TeV is possible with 300 fb -1 , and as little as 1 fb -1 is needed to probe 
m H ±± < 400 GeV. 

Discovery of H ±:k with m#±± < 400 GeV would enable precise measurements of the branch- 
ing ratios (BRs) of H ±:iz — > if if with the anticipated final integrated luminosity at the LHC. 
Models which predict BR(if ±± — > if if) in terms of experimentally constrained and/or mea- 
sured parameters are of particular phenomenological interest. In the Higgs Triplet Model 
(HTM) [3], [4] neutrinos acquire a Majorana mass given by the product of a triplet Yukawa 
coupling (hij) and a triplet vacuum expectation value va- Consequently in the HTM there 
is a direct connection between and the neutrino mass matrix which gives rise to phe- 
nomenological predictions for processes which depend on h^. Since the coupling determines 
BR(i7 ±=t — > if if), this mechanism of neutrino mass generation can be tested if precise measure- 
ments of BR(iJ ±=t — » if if) are available [6,7]. A detailed quantitative study of the dependence 
of hij on all the neutrino oscillation parameters has not yet been performed (for previous anal- 
yses see [7, 17]). 

Of particular interest is the dependence of on the absolute neutrino mass and Majorana 
phases which is the focus of the present work. Those parameters, which cannot be probed in 
neutrino oscillation experiments, would significantly affect BR(iJ ±± — > if if). We perform a 
study of the capability of the LHC to probe the neutrino mass spectrum and Majorana phases 
assuming that neutrino mass is generated solely by the combination h^v^. in the HTM. In 
particular, we investigate the possibility to establish m ^ and/or CP- violation from Majo- 
rana phases at the LHC by means of a x 2 analysis with three 21 channels of if ±± decays. It is 
extremely difficult for neutrinoless double beta decay experiments [18] to measure CP-violation 
from Majorana phases because they affect this process in combination with unmeasured pa- 
rameters [19], while the absolute neutrino mass can only be measured directly by the future 
Tritium beta decay experiment [20] if m > 0.2eV. 

Our work is organized as follows: in section 2 we briefly review the HTM; section 3 describes 
the phenomenology of H =t± at hadron colliders; the numerical analysis is contained in section 
4 with details of a x 2 analysis presented in the appendix; conclusions are given in section 5. 
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2 The Higgs Triplet Model 



In the Higgs Triplet Model (HTM) [3], [4] a 7 = 1, Y = 2 complex SU(2) L isospin triplet of 
scalar fields is added to the SM Lagrangian. Such a model can provide a Majorana mass for the 
observed neutrinos without the introduction of a right-handed neutrino via the gauge invariant 
Yukawa interaction: 

C = hij^f L CiT 2 Aif}j L + h.c (1) 

Here hij(i,j = 1, 2, 3) is a complex and symmetric coupling, C is the Dirac charge conjugation 
operator, t 2 is a Pauli matrix, ipn = (V«, li) T L is a left-handed lepton doublet, and A is a 2 x 2 
representation of the Y = 2 complex triplet fields: 



A° -A+/V2 

A non-zero triplet vacuum expectation value (A ) gives rise to the following mass matrix for 
neutrinos: 

rriij = 2hij(A°) = V2hijV A (3) 

The necessary non-zero va arises from the minimization of the most general 577(2) (g> U(1)y 
invariant Higgs potential, which is written as follows [6,7] (with $ = (0 + ,0°) T ): 

V = m 2 ($ f $) + Ai($ t $) 2 + M 2 Tr(A t A) + A 2 [Tr(A f A)] 2 + A 3 Det(A t A) 



+A 4 ($ t $)Tr(A t A) + A 5 ($ t r i <l>)Tr(A t r i A) + \-^=p{® T ir 2 A^) + h.c^j 



(4) 



Here m 2 < in order to ensure (0°) = v/y/2 which spontaneously breaks SU(2) (g> U(l)y 
to U(1) Q , and M 2 (> 0) is the mass term for the triplet scalars. In the model of Gelmini- 
Roncadelli [21] the term p{^ T ir 2 A^) is absent, which leads to spontaneous violation of lepton 
number for M 2 < 0. The resulting Higgs spectrum contains a massless triplet scalar (majoron, 
J) and another light scalar (H°). Pair production via e + e~ — > H°J would give a large contri- 
bution to the invisible width of the Z and this model was excluded at LEP. The inclusion of 
the term /i($ T ^A^Q) explicitly breaks lepton number when A is assigned 7 = 2, and elimi- 
nates the majoron [3,4]. Thus the scalar potential in eq. (prj together with the triplet Yukawa 
interaction of eq. (PQ) lead to a phenomenologically viable model of neutrino mass generation. 
The expression for v& resulting from the minimization of V is: 

v A ~ fxv 2 /2M 2 (5) 

In the scenario of light triplet scalars (M w v) within the discovery reach of the LHC, eq. (151) 
leads to va ~ A 4 - In extensions of the HTM the term /x($ T ir2A^$) may arise in various ways: 
i) the vev of a Higgs singlet field [22,23]; ii) be generated at higher orders in perturbation 
theory [7]; iii) be generated in the effective Lagrangian [10]; iv) originate in the context of extra 
dimensions [6,8]. 

An upper limit on va can be obtained from considering its effect on the parameter p{= 
Mw/Mz cos 2 9w)- In the SM p = 1 at tree-level, while in the HTM one has (where x = va/v): 

1 r l + 2x 2 
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The measurement p « 1 leads to the bound v&/v <^0.03, or v a < 8 GeV. At the 1-loop level 
v a must be renormalized and explicit analyses lead to bounds on its magnitude similar to those 
derived from the tree-level analysis [24]. 

The HTM has seven Higgs bosons (H ++ , H , H + , H~, H°, A , h°). The doubly charged 
if ±± is entirely composed of the triplet scalar field A ±=t , while the remaining eigenstates are 
in general mixtures of the doublet and triplet fields. Such mixing is proportional to the triplet 
vev, and hence small even ifv& assumes its largest value of a few GeV. Therefore H ± ,H°,A° 
are predominantly composed of the triplet fields, while h° is predominantly composed of the 
doublet field and plays the role of the SM Higgs boson. The masses of H ±± ,H ± ,H°,A° are 
of order M with splittings of order \$v . For M < 1 TeV of interest for direct searches for the 
Higgs bosons at the LHC, the couplings hij are constrained to be 0(0.1) or less by a variety of 
processes such as \i — > eee, r — > III etc. which are reviewed in [25,26]. 



3 Production and Decay of H ±zL at Hadron Colliders 

The most distinct and experimentally accessible decay mode of if ±± is to two same-sign charged 
leptons [27]. Without loss of generality one can work in the basis in which the charged lepton 
mass matrix is diagonal i.e., if are the mass eigenstates. Then the decay rate for BR(if ±=t — > 

l t l t) (h3 = e > T ) is § iven b y : 



r(H±± it$) = (7) 



where S — 1(2) for i — j (i j). Clearly T(H ±i: — > if if) depends crucially on the absolute 
values of the h^, where is related to the neutrino mass matrix via eq. (j3J). However, if no 
other decay modes for if ±± are open kinematically the leptonic BRs are determined solely by 
the relative values of h^. Other decay modes for H ±zL can be important in regions of parameter 
space e.g., i) if ±=t — > H ± W*, which is potentially sizeable for m H ±± > m H ± [17,28], and ii) 
if ±± — > W ± W ± , which is proportional to the triplet vev, In this work, we assume > 
rriH±± (which precludes H ±zL H^W*) and v a < 1 MeV, which suppresses if ±± — > W^W^ 
sufficiently in the HTM (e.g., see [29]). Then, -»• if if can be regarded as the sole decay 
mode for if ±± and one has: 



BR W / = BR(iy ±=b — >• if if) = v^nr^fi- (S) 

ij I'M 



3.1 Searches for i7 ±=b at the Tevatron 

In the year 2003 the Fermilab Tevatron performed the first search for H ±zSz at a hadron collider ID 
The DO collaboration [13] searched for if ±± — > fi // while the CDF collaboration [14] searched 
for 3 final states: if ±± — > e ± e ± , e ± /i ± , /j, =fc /j- ± . The assumed production mechanism for i7 ±± is 
qq — > 7*,Z* — > H ++ H , which proceeds via gauge strength couplings and depends on only 
one unknown parameter, m#±± [25,32]. 



Direct searches for H ±± have also been performed at LEP [30] and HERA [31]. 
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e/x 


/i/i 


er 


fiT 


TT 


21 


> 133 GeV 


> 113 GeV 


> 136 GeV 








31 








> 114 GeV 


> 112 GeV 




41 








> 114 GeV 


> 112 GeV 





Table 1: Mass limits on m#±± from searches for H ±zSz — > Z^Z^ at Tevatron Run II. 




The searches performed in [13,14] seek at least one pair of same-sign leptons with high invariant 
mass i.e., the search is sensitive to single production of if ±=t . The SM background can be reduced 
to negligible proportions with suitable cuts. Single if ±± production mechanisms which involve 
a dependence on potentially small parameters such as the Yukawa coupling hij or triplet vev v& 
are subdominant at Tevatron energies (e.g., see [33]). In [17] it was suggested that this search 
strategy is also sensitive to the mechanism q'q — > W* — > if ±=t if T [34] , which has a cross-section 
comparable in magnitude to that of qq — > 7* , Z* — > H ++ H . The following inclusive single 
if ±± cross-section was introduced, which would extend the search sensitivity to larger values 
of rriH±± and strengthen the mass limits on m#±± derived in [13, 14]: 

a H ±± = a(qq -> 7*, Z* -> ) + 2a(gg -> -> H ++ H') (9) 

Here the factor of 2 accounts for the CP conjugate process — > — >■ H H + . 

In 2006 the CDF collaboration searched for if ±=t decays involving t^ [35]. The strategy of 
searching for one pair of same-sign leptons (2/) is not effective due to the larger SM backgrounds, 
and instead three (3/) and four (4/) lepton searches were performed. The production mechanism 
qq — > 7*, Z* — > H ++ H was assumed. The process qq — > VT* — > H ±:iz H T never contributes to 
the 4/ signature, but can contribute to the 3/ signature if H ± decays leptonically, H ± — > l^v. 

In Table [TJ the mass limits for m#±± from the Tevatron searches are summarized. A blank 
entry signifies that no search has yet been performed. The displayed mass limits assume 
production via qq — > 7*, Z* — > H ++ H for a H ±:k belonging to a SU(2)l triplet with Y = 2. 
Moreover, BR( J ff ±± -> Zf Z±) = 100% in a given channel is assumed. The ultimate sensitivity 
at the Tevatron is expected to be m^±± ~ 250 GeV in the ee, e\x and \i\x channels. 

3.2 Simulations of H ±zt production at the LHC 

Several simulations have been performed for if ±± — > ifl^ = e, /x, r) at the LHC [15,16, 
29,36,37]. The production mechanism is assumed to be qq — > 7*,Z* — > H ++ H followed 
by H ++ H — > llll. The LHC sensitivity to — > ifl^ considerably extends that at the 
Tevatron due to the increased cross-sections and larger luminosities e.g., the analysis of [16] 
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m H ±± (GeV) 


N Al (30 fb" 1 ) 


N 4l (300 fb" 1 ) 


N 2l (300 fb" 1 ) 


200 


1500 


15000 


42000 


300 


300 


3000 


8400 


400 


90 


900 


2500 



Table 2: Approximate number of events for pair production of ff ±± (N41) and single production 
of if ±=t (N 2 i) at the LHC for integrated luminosities of C = 30 fb -1 and C = 300 fb -1 with 
efficiency e c s = 0.5. We assumed m H ±± = mn± to calculate N 2 i- 

shows that a if ±=t can be discovered for m H ±± < 800 GeV assuming BR(if ±± — > fi^fi^) = 100% 
and £ = 50 fb -1 . Importantly, all the above simulations suggest that as little as 1 fb -1 is needed 
for discovery of m H ±± < 400 GeV if one of BR(iJ ±± — > e ± e ± , e^fi^, /i ± /i ± ) is large, and hence 
such a light if ±± would be found very quickly at the LHC. 

The sensitivity of the LHC to single production of if ±± — > if if for i,j = e, fx has only 
been performed in [37], and importantly the SM background was shown to be negligible in 
the signal region of high invariant mass. It was concluded that such a search strategy allows 
more if ±± — > if if events than the 4 lepton search since the event number is linear (and not 
quadratic) in BR(if ±± — > if if)- Therefore the 21 search is more effective at probing small 
BR( J ff ±± -> If if). Importantly, the addition of the channel qq -> W* -> H^H^ (eq. (JSJ)) 
would further enhance the event number for a given m#±± . 

In Table |2] we show approximate expected numbers of 21 and 41 events arising from pair 
and singly produced H ±zSz — > if if at the LHC. We only consider the decay channels H ±zSz — > 
e ± e ± , e ± /i ± , which offer the greatest if ±=t discovery potential. A detection efficiency of 

0.5 is assumed, which is slightly less than the values given in [16] for if ±=t — > /i ± /x ± . The 
theoretical if ±=t cross-section is multiplied by this detection efficiency and the SM background 
is taken to be negligible. The number of 4/ events for a specific m#±± is denoted by N41, 
assuming integrated luminosities of C = 30 fb -1 and C = 300 fb -1 . The displayed numbers are 
for BR(ii ±± -> Zf = 100% in a given channel. In the HTM, BR(ii ±± -> if if) is necessarily 
< 100% (eq. © and eq. ©) and hence N 4l must be multiplied by [BR( J ff ±± -> Z^Z^)] 2 . The 
final column shows the number of 2Z events (A^z) obtained by adding the contribution from the 
mechanism qq — > W* — > iJ ±± if =F as defined in eq. ([9]). We take m H ±± = m H ± which increases 
the number of singly produced if ±± events by a factor of around 2.8 for 200 GeV < m#±± < 
400 GeV [17]. For BR( J ff ±± -> if if) < 100%, the numbers presented in Table [2] are scaled as 
shown in Appendix lAl 

It is clear from Table [2]that early discovery of if ±± at the LHC with ttih±± < 400 GeV would 
allow large event numbers for if ±± with the expected integrated luminosities of £ = 300 fb -1 . 
This would enable precise measurements of BR(if ±=t — > e ± e ± , e ± /i ± , /x^// 1 * 1 ) for the dominant 
channels. Sensitivity to BR(iJ ±=t — > e ± e ± , e ± /i ± , ^fi^ ~ 1% or less would also be possible in 
the 2Z channel. 
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4 Numerical Analysis 



The mass matrix for three Dirac neutrinos is diagonalized by the MNS (Maki-Nakagawa-Sakata) 
matrix Vmns [38] for which the standard parametrization is: 

C12C13 S12C13 s 13 e~ lS 

V MNS = I - S 12C23 - Ci 2 S23Si 3 e l<5 Ci 2 C 23 ~ Si 2 S23Si 3 e i<5 S 23 Ci 3 | , (10) 

S12S23 ~ c 12 c 2 3Si 3 e lS -c i2 s 23 - si 2 c 23 si 3 e Jl5 c 23 ci 3 

where s»j = sin 6^- and Cy = cos 6^, and 5 is the Dirac phase. For Majorana neutrinos, two 
additional phases appear and then the mixing matrix V becomes 

V = Vmns x diag(l, e^' 2 , e^ 2 ), (11) 

where ipi and tp 2 are referred to as the Majorana phases [3,39]. Since we are working in the 
basis in which the charged lepton mass matrix is diagonal, the neutrino mass matrix is then 
diagonalized by V. Using eq. fl3j) one can write the couplings hy as follows [6,7]: 

hij = [VMNsdiag(mi,m 2 e il/ ' 1 ,m 3 e 4lp2 )VMNs]i i ( 12 ) 



Then, eq. (jH) becomes 



S { [VWsdiagK, m 2 e^, m 3 e i ^)V^ NS ] j 
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BR(H -+ Iflf) = =—5 (13) 



Note that the branching ratios are independent of m^±± and v a, and given by neutrino pa- 
rameters only. 

Neutrino oscillation experiments involving solar [40], atmospheric [41], accelerator [42], and 
reactor neutrinos [43] are sensitive to the mass-squared differences and the mixing angles, and 
give the following preferred values: 

Am^ = m\ - ml ~ 7.9 x 10 _5 eV 2 , \Am 2 31 \ = \m\ - m\\ ~ 2.7 x 10" 3 eV 2 , (14) 
sin 2 26 12 ~ 0.86 , sin 2 2^ 23 ~ 1 , sin 2 29 13 ^ 0.13 . (15) 

The small mixing angle #i 3 has not been measured yet and hence the value of 6 in completely 
unknown. Since the sign of Am^ is also undetermined at present, distinct neutrino mass 
hierarchy patterns are possible. The case with Am 2 x > is referred to as Normal hierarchy 
(NH) where m\ < m 2 < m 3 and the case with Am 2 x < is known as Inverted hierarchy (IH) 
where m 3 < m\ < m 2 . Information on the mass mo of the lightest neutrino and the Majorana 
phases cannot be obtained from neutrino oscillation experiments. This is because the oscillation 
probabilities are independent of these parameters, not only in vacuum but also in matter. If 
m > 0.2eV, future 3 H beta decay experiment [20] can measure it. Experiments which seek 
neutrinoless double beta decay [18] are sensitive to only a combination of neutrino masses and 
phases. Certainly, extracting information on Majorana phases alone from these experiments 
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seems extremely difficult, if not impossible [19]. Therefore it is worthwhile to consider other 
possibilities. 

Multiplying out eq. (|T2|) one obtains the following explicit expressions for hif 

U 1 ( 2 2, 22 ii£i, 2 -2i<5 itpo\ 

h ee = ^ ^iCi2 c i3 + m 2S 12 c 13 e Vl + m 3 s l3 e e V2 j , 

Kf, = |mi(-Si 2 C23 - Ci 2 g23-Sl3e t<5 )ci2Ci3 

+ rn 2 (ci 2 c 23 - si 2 S23Si3e J,5 )si2Ci3e iipi + m 3 s 23 c 13 si 3 e~ lS e ttp2 1 , 

ft-er — | r ^l('Sl2'S23 — c 12C23'5l3 e *' 5 ) c 12Cl3 

+ ^2(-ci2S 2 3 - s 12 c 23 s 13 e lS )s 12 c 13 e llfl + m 3 c 23 Cx 3 Si 3 e~ lS e lt(>2 } , 
h MM = — |mi(-si 2 c 2 3 - ci2S23Si3e* 5 ) 2 + m 2 (ci 2 C23 - s 12 s 23 s 13 e t5 ) 2 e t<Pl + m 3 s 2 23 c\ 3 e lip2 ^ , 
hpr = \rn x {-s X2 c 23 - c 12 s 23 s 13 e lS ) (S12S23 - c 12 c 23 s 13 e lS ) 

+ m 2 (c 12 c 23 - si 2 s 2 3S 13 e lS )(-c 12 s 2 3 - snc^s^e^e^ 1 + m 3 c 23 s 23 c\ 3 e llf2 ^ , 

Kt = — |^i(si 2 S23 - c 12 c 2 3S 13 e tS ) 2 + m 2 (-c 12 s 23 - si 2 c 2 3Si3e l& ) 2 e llfl + m 3 cl 3 cl 3 e t<P2 ^ . 

(16) 

One may express mi,m 2 ,m 3 in terms of two neutrino mass-squared differences (Am^Am^) 
and the mass of the lightest neutrino mo- The hij are functions of nine parameters: 
Am^Am^mo, three mixing angles (#12, #13, #23) and three complex phases (5, ipi, <p 2 ). Four of 
the above parameters have been measured well experimentally (see eq. IT5I) and thus the HTM 
already provides numerical predictions for BR(iJ ±=t — > if if) as a function of the unmeasured 
five parameters. Four cases corresponding to no CP violation from Majorana phases can be 
defined as follows: Case I (tp\ — 0, <p 2 — 0); Case II (<pi = 0, ip 2 = 7r); Case III (<pi = tc, tp 2 = 0); 
Case IV ((fx — ir, <p 2 — tt) • These four cases have been studied in [7, 17] for values of m = 
or O(l) eV. In this work we will quantify the dependence of on m , (p% and <p 2 . Those 
3 parameters essentially determine BR(if ±=t — > if if), with subdominant corrections from the 
neutrino oscillation parameters. Hence multiple signals for H ±zSz — > if if at the LHC would 
probe mo, tpi and ^2, in the context of the HTM. 

We note here that such collider probes of the Majorana phases are particular to models in 
which lepton number violation (which leads to the Majorana neutrino mass) is associated with 
New Physics particles at the TeV scale. Analogous collider probes are not possible in models 
where the scale of lepton number violation is much higher e.g., supersymmetric models with 
very heavy right-handed neutrinos with masses of order 10 12 GeV. However, in such models the 
Majorana phases (which are also required for successful leptogenesis) can significantly affect 
the rates for the lepton flavour violating (LFV) decays r — ► Z7 and fi — > e'-f [44]. Likewise, in 
the HTM ipi and ip 2 would also affect the rates for LFV decays which depend explicitly on 
e.g., r — > III and \i — > [7]. Similar studies of the effects of Majorana phases on these LFV 
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Figure 1: BR(if ±± — > if if) with no CP violation from Majorana phases (ipi = <^2 = 0) as a 
function of mo in NH (left) and IH (right). 



decays in other models which contain a TeV scale H and an analogous coupling have 
been performed in [45]. 

4.1 Dependence of BR(i/ ±± -> / ± / ± ) on the neutrino mass spectrum 

In Fig. [H we plot BR(if ±± — > Z 4 L ) as a function of m for cpx — and <£> 2 = (Case I) in NH 
and IH. The values of the oscillation parameters are fixed in the figures as follows: 

Am 2 21 = 7.9 x l(T 5 eV 2 , \Am 2 31 \ = 2.7 x 10~ 3 eV 2 , sin 2 29 12 = 0.86 , (17) 
sin 2 2fl 23 = l, sin 2 2^13 = 0.13, 6 = 0. (18) 

For NH with smaller m , the e-related modes are suppressed because mi{= m ) and m 2 
are small and the contribution from the heaviest mass m 3 is also small due to the tiny 6 13 . 
BR^, BR MT , and BR TT are roughly equal El, and they dominate for smaller m because 
does not appear with #i 3 . Note that BR M/J ~ BR TT and BR eAt ~ BR er can be understood as the 
approximate symmetry for /x-r exchange by virtue of sin 2 2# 2 3 = 1 and tiny # 13 . For larger mo, 
BR(i7 ±=l= — > if If) dominate because all coefficients of m; are positive for ha with = in 
Case I and there is no strong cancellation among these terms even for #i 3 ^ 0. For mo ^ 0.3 eV, 
the results for NH and IH are almost identical because of almost degenerate masses. Since BR 
are given by ratios of \hij\ 2 , they converge for large mo- This is an attractive feature because 
HTM can predict certain ranges of BR without restricting mo- 

On the other hand, IH case gives rather simple results. This is because mi and m 2 include 
the larger scale yjAm|J without any suppression from 13 and consequently the effects of small 
Am 21 and 813 are hidden. In this case, BR ee becomes the dominant channel even for smaller 
mo; one naively expects BR ee ~ 4BR WI for m = with eq. (|T6l) . 

2 The naive expectation is that BR MT is twice as large as the other two modes because of |^t| 2 + |ft-r^| 2 - 
However, this difference is accidentally compensated by the effect of Am^. 
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Figure 2: BR(iy ±=l= — > e ± e ± ) as a function of mo in NH (left) and IH (right). We take four 
cases of CP conservation (I-IV) for sin 2 29r 3 = 0.13 and 4 values of 5 (= 0, 7r/2, tt, 3tt/2). The 
dotted lines are obtained for #13 = 0. 



In order to quantify the effect of the unmeasured (VjviNs)e3 on BR/;/, we show BR ee for 
sin 2 2^13 = 0.13,0, and 4 values of the Dirac phase 5 (= 0, tt/2, it, 3ti/2) for Cases I to IV in 
Fig. [2j Other parameters are same as Fig. [TJ Several lines are coincident in the figure, since 
5 and ifi2 always appear as a combination 25 — ip2 in h ee . For example, Case I with 5 = 
and Case II with 5 = n/2 give the same lines. Although the contribution of (T4iNs)e3 to BR ee 
is considerably smaller than the effect of varying mo or the Majorana phases (as expected by 
quadratic suppression with small 9 13 ), it is not negligible. Consequently, it is also not negligible 
for other the BR. Fig. [2] also shows that the HTM predicts small BR ee in case III and IV (cases 
with if \ = 7r) for any value of mo, #13, and 5 in both of NH and IH. Thus, Fig. [2] indicates 
that some information on Majorana phases may be extracted without knowledge of the sign of 
Amjj and the values of m , #13, and 5. 

4.2 Dependence of BR(i/ ±± -> if if) on Maj orana phases 

In this section we show the dependence of BR(if ±± — > if if) on Majorana phases with the 
values of the oscillation parameters given in (1171) and (ITS]) . For m = 0, only the relative phase 
<A-ei(= (pi — P2) determines in NH. In Fig. [3] (left) we plot BR(if ±± — > iflf) as a function 
of ip ie \ for NH with m = 0. Tiny 9i 3 suppresses the dependence of BR ee on <p Te \, and then 
the HTM gives a clear prediction of very small BR ee for this case. The other BKw change 
non-negligibly e.g., BR eA1 and BR W vary by 0.05 or more, which could be larger than the 
experimental error if sufficiently large numbers of H ±zL are produced. In IH for mo = 0, BR/// 
does not depend on ip 2 because it always appears multiplied by m . Fig. [3] (right) shows the 
dependence of BR/// on tp\ for this case. One can see that the dependence on (fx is even more 
pronounced because m 2 for m = in IH (= yfAm|J) is larger than that in NH (= a/ Am^). 
In particular, BR ee and BR eM seem to be very useful for extracting information on ipx because 
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they have large and opposite dependence. This dependence on <pi can be understood by the 
relative sign of the terms of mi and 777.2 in eq. (fTBT) . which is + for h ee and — for h efl , neglecting 

#13- 

For mo = 0.3 eV, where neutrino masses are almost degenerate, Fig. H] shows that the phase 
cpi has a large effect on all BR(if ±i: — > if if). The figure also shows that the (f2 dependence 
for BR eM and BR er is sizeable and overcomes the suppression from #13. However, the latter 
suppression factor ensures that BR ee has a very small dependence on <p 2 , and so this channel is 
crucial for extracting clear information on tp\ and/or m without contamination from the effect 
of y? 2 - 

It is clear that the Majorana phases have a large effect on BR/// in the HTM, and their 
inclusion is required in order to quantify the allowed regions of BR///. 



4.3 Sensitivity to sigr^Am^) and ijiq. 

As we have seen, both the neutrino mass spectrum and the Majorana phases have large effects on 
BR(i/ ±± -> if if) in the HTM, but the model also gives a clear prediction for BR/// . Therefore 
we can expect to obtain some information on the neutrino mass spectrum and/or the Majorana 
phases by observing BR///. In our analysis, we use BR(if ±± — ► e ± e ± ,e ± fi ± , ^fi^) for which 
the LHC expects greatest sensitivity. Naively, three measurements of BR/// are sufficient to 
extract information on the three parameters. 

In this section we consider the possibility to determine sign(Am| 1 ) and/or to exclude tuq = 
0. Hereafter, we use (fT71) and 

sin 2 29 23 > 0.93 , sin 2 29 13 < 0.13 , 5 = 0-2n. (19) 

Non-zero values of #13 and 5 affect BR ee as was shown in Fig. [2j Moreover, deviation of #23 from 
7r/4 especially affects BR^ in our analysis because a rather wide range 0.37 < S23 < 0.63 is 
allowed. In Fig. the allowed regions of BR/// in the HTM are shown by the shaded (light and 
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NH-IH 



NH-IH 




dark) regions; the dark shaded area corresponds to the overlap of the allowed regions for NH 
and IH. The area above the dotted line is unphysical because the sum of BR exceeds unity. It 
is clear that the HTM predicts BR ee < 0.49, BR efM < 0.61, and BR wt < 0.47. The areas outside 
the solid line are unreachable, and are particular to the HTM (in which neutrino mass is solely 
given by eq. ([3])) and will differ from the corresponding disallowed regions, if anyjf] in other 
models which can contain a light if ±± . Hence in the event of signals for if ±± — > if if these 
allowed regions can be compared with the experimental measurements of BR(if ±=t — > if if) in 
order to test the HTM. If measured BR are outside the shaded region the HTM is disfavored. 
If given sie nals for -> if if are compatible with the reachable regions in the HTM, the 
data can be interpreted in the context of this model and one can try to extract information on 



3 If hij are arbitrary parameters, as in the Left- Right symmetric model of [1], the unreachable region would 
vanish. 
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Figure 5: The shaded area surrounded by the solid line is attainable in the HTM. Allowed 
regions for NH and IH correspond to BR inside of the dashed and dash-dotted lines, respectively. 
The dotted line divides the physical and unphysical regions. 



parameters of neutrino mass matrix. 

Furthermore, the area inside the dashed and dash-dotted lines in Fig. [5] is possible for NH 
and IH, respectively. NH gives the additional bounds BR ee < 0.33 and BR W < 0.51, while IH 
gives 0.03 < BR ee and BR wt < 0.34. These are also particular features of the HTM. If BR in 
the light shaded region are measured, NH or IH is disfavored in the HTM. 

Next, let us consider the possibility to exclude m = 0. Fig. [H] shows the attainable regions 
of BR with m = 0. As in Fig. the area inside the thin solid line is reachable in the HTM, 
and the region above the dotted line is unphysical. The dashed line corresponds to the BR 
which can be obtained in NH, and the dash-dotted line is BR in IH, both with mo = 0. Note 
that the area inside these lines can also be obtained with m 7^ but outside is impossible for 
m = 0. This behaviour can also be seen in Fig. [3J The bold solid line denotes a x 2 analysis for 
measurements of signals in the 21 channels (N ee , N efl , and N^) of if ±=t decays (see Appendix IA1 
for detail). If BR in the shaded region are measured, one can exclude mo = at 90%CL in 
HTM. NH with mo = gives a very clear prediction of small BR C ™ as can be seen in Fig. [3j 
The bold line shows that mo = can be excluded at 90%CL if BR*™ 6 lies in the narrow interval 
0.02 < BR*™ 6 < 0.05, even without measuring other BR. If only BR efl or BR MM is observed it 
is impossible to exclude m = 0. Therefore, it is essential to use more than two BR to exclude 
m = 0. 



4.4 Sensitivity to CP violation from Majorana phases 

In this section, we discuss the possibility to exclude the CP-conserving cases for Majorana 
phases, namely <£>i,<£>2 = 0, 7r. Establishing experimentally that CP is violated by Majorana 
phases in the HTM is of much phenomenological interest and is feasible if there are sufficiently 
large numbers of if ±=t . 

In Fig. the attainable regions of BR for the CP-conserving cases of Majorana phases 
are shown. As in Fig. [SJ the area inside the thin solid line is reachable in the HTM, and the 
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Figure 6: The thin solid line shows the attainable region in the HTM and the dotted line 
divides the physical and unphysical regions. The area inside the dashed and dash-dotted lines 
can be obtained with mo = in NH and IH, respectively. The bold solid line is a result of a x 2 
analysis and m = can be excluded if BR are measured within the shaded region. 



region above the dotted line is unphysical. Bold dashed, thin dashed, bold dash-dotted, and 
thin dash-dotted lines show the possible regions for Case I, II, II, and IV, respectively. From 
the figures it is clear that the HTM predicts very specific regions for BR( J ff ±± If If) in the 
CP-conserving cases, especially in the space of BR ee -BR eA1 ; BR ee is small (< 0.07) in Case III 
and IV (cases with <p\ = n) as seen in Fig. [2J and BR e/i cannot be so large (< 0.09) in Case I 
and II (cases with (pi = 0). The bold solid line is obtained by a \ 2 analysis with 21 channels of 
H ±i: decays (see Appendix [A] for details), and the CP-conserving cases are excluded at 90%CL 
in the HTM if experiment measures BR in the shaded region. 

4.5 Accuracy of the determination of uiq and Majorana phases 

Once measurements of BR ee , BR eM , and BR W are available one can constrain mo and the 
Majorana phases. In this section we show how precisely they can be constrained by assuming 
BR ee = 0.268, BR eM = 0.243, and BR M/J = 0.066 as an example of a possible signal. These 
values can be obtained by taking m — 0, <pi — ir/2, sin 2 26 ) 13 = 0.1, 5 = 0, and # 2 3 = 7r/4 for 
IH in HTM. It is evident from Fig. that these values of BR can be accommodated in the 
HTM, but not in the case of NH - see Fig. [5](a). 

In Fig. [HI the allowed regions at 90%CL are shown in the spaces of v?i-y?2, rno-fi, an d ^o-y?2- 
They are obtained by a \ 2 analysis with three numbers of 21 channel signals (N ee , N e/1 , and 
iV^) of H ±i: decays (See Appendix [A] for details). Since the above BR are inside the shaded 
regions in Fig. [71 Case I-IV of CP-conservation are all excluded which is clearly depicted in 
Fig. Eta). Two regions in Fig. [8](a) are just copys of each other because BRw are unchanged by 
replacing hij with h*j. Although (fx can be constrained very well, all values of ip 2 are allowed 
in Fig. [8] because BR in IH with mo = does not depend on (p 2 . In general one expects looser 
constraints on cp 2 because dependence of BR ee and BR eA1 on (p 2 is suppressed by tiny 9 ri . Fig. [S] 
also shows that a stringent upper bound on m Q (m < 0.07 in this example) can be obtained 
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Figure 7: The thin solid line shows the attainable region in the HTM and the dotted line 
divides the physical and unphysical regions. Cases I, II, III, and IV with CP conservation 
for Majorana phases can give BR within the bold dashed, thin dashed, bold dash-dotted, and 
thin dash-dotted lines, respectively. The bold solid line is a result of a y 2 analysis and the 
CP-conserving case is excluded at 90%CL if BR within the shaded region are measured. 



by measuring BR. 



5 Conclusions 

Doubly charged Higgs bosons (i? ±± ) are a distinctive signature of a variety of models which 
can accommodate neutrino mass. Discovery of if ±± at the Tevatron/LHC with a reasonably 
light mass (< 400 GeV) could enable thousands of H ±±L — > l { I- events to be observed with the 
anticipated final integrated luminosities at the LHC. In such a scenario precise measurements 
of if ±=l= — > if if would be feasible and would constitute a means of experimentally favouring or 
disfavouring various candidate models which contain a H ±:k . 

The Higgs Triplet Model (HTM), in which (Majorana) neutrino mass is solely given as 
the product of a triplet vacuum expectation value va and a triplet Yukawa coupling hij, is a 
particularly simple model which leads to a specific relationship between the branching ratios 
(BR) of H ±:k -> if if and the neutrino mass matrix. In particular, BRs are determined only 
by the neutrino mass matrix independently of va and Higgs boson masses if m#±± < m#±. 
We performed a quantitative study of the dependence of the BRs on all the parameters in the 
neutrino mass matrix, especially the dependence on the three parameters which induce the 
most uncertainty in such predictions: the lightest neutrino mass m and the two Majorana 
phases. We displayed the allowed regions of BR(ii" ±± — > if if) for i, j = e, fx and showed that 
they are considerably smaller than for the case of an arbitrary . By measuring these BRs it 
is possible to both test the HTM and determine the sign of Am^. 

Due to the large effect of mo and the Majorana phases on the BRs we showed that it 
is possible to extract information on these three parameters by measuring BR ee , BR eAt , and 
BR MM . Such information cannot be obtained from neutrino oscillation experiments. By using 
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Figure 8: The shaded region is allowed at 90%CL for BR ee = 0.268, BR eM = 0.243, and 
BR^ = 0.066. These values can be obtained by taking m = 0, ipi = tt/2, sin 2 26> 13 = 0.1, 
S = 0, and 6 23 = vr/4 for IH in HTM. 



the 21 channel of H ±i: decays, it was found that m = and/or the CP-conserving case for 
the Majorana phases can be excluded without knowledge of sign(Am| 1 ), m , #13, and 5. There 
is a small range of BR ee (0.02 < BR ee < 0.05) where m = can be excluded at 90%CL 
by measuring only BR ee . For the exclusion of the CP-conserving case for Majorana phases, 
BR ee and BR eM are especially important. For a specific choice of BR ee , BR eM , and BR W we 
displayed the allowed regions of (m , ipi, ^2) and showed that it is possible to constrain m and 
(pi stringently. 

In contrast, establishing that Majorana phases are non-zero is notoriously difficult in models 
in which lepton number violation is associated with particles at a very high energy scale (e.g. 
right-handed neutrinos with a Majorana mass of order 10 10 GeV or more). In the HTM with 
a light if ±± this becomes is a realistic and unique possibility. 

Note added: During finalization of this paper an article [46] appeared which deals with the 
same topic. Although there is inevitably some overlap, our analysis differs considerably from 
theirs e.g., i) our definition of the x 2 is a function of three BRs in the 2 lepton channel, while 
their \ 2 is f° r fi ye BRs in the 4 lepton channel, and ii) the presentation of results, especially 
figures 5-8. Where comparisons can be made we find good agreement with their results. 
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A Definition of A% 2 

In order to test assumptions, e.g., CP-conservation of Majorana phases, we calculate a Ax 2 . 
In our analysis, BR*™ 6 , BR*" 16 , and BR*™ 6 are used as input which can be understood as true 
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values chosen by nature or best fit values from experiment. We try to fit these true values with 
BRf^ which are calculated theoretically by eq. (|HJ) for the assumption to be tested. We utilize 
the 21 channel and the number of events is given by 

N u > = e eff {iV pair BR«'(2 - BR W ) + iV singlc BR^} , (20) 

where iVp air (iV S i n gi e ) is the number of pair (single) production of doubly charged Higgs, and e e s 
denotes an efficiency due to event cuts; we can obtain the numbers of i7 ±± production from 
Table [2] as iV pair = N 4l and iV sing i e = N 2[ - N it . We take e efr = 0.5. We use iV pair = 900 and 
-^single = 1600 by assuming m H ±± = m H ± for simplicity (see Table [2]) and because the difference 
of their masses should not be too large in order to maintain p ~ 1 at the 1-loop level. For the 
number of remaining background events after cuts, we use Nbg — 1 f° r each Nu> because the 
signal is expected to be almost background free. In the fitting procedure, A% 2 is minimized 
within possible regions of the following parameters: 

m f = 0-0.3eV, iff ,(pf = 0-2vr, (21) 
sin 2 20** = 0-0.13, 5 fit = 0-2tt, (22) 
sin 2 2dl\ > °- 93 , hierarchy 61 = NH, IH. (23) 

The minimization gives the most pessimistic value for exclusion of the assumption because a 
smaller value of Ay 2 means a better fitting with the assumption. The other parameters are fixed 
for true and fitting values as ffTTj) . Since we wish to extract information on three parameters 
(mo, (fx, (f2), our analysis is for three degrees of freedom and Ay 2 = 6.3 corresponds to 90%CL 
even if results are projected onto spaces of two parameters. 



A.l Exclusion of mo = and CP-conserving case for Majorana phases 

In Fig. [6l the region where mo = is excluded in the HTM is given by calculating the Ax 
function. The Ay 2 is defined as 



AxL(BR t e r e ue ,BR t e r u ue ,BR t u r u ue ) = min I V \f -p- 



(24) 



x = {On, 6, 9 23 , hierarchy} . (25) 

where A^ uc is given by (BR*™ , BR*™ e , BR*™ 6 ) and Np is calculated with the theoretical form 
of BR (eq. (jSJ)) for {i^,i^,m^,x ). When we show the projected region of Ax 2 mo < 6.3 
onto the plane BR*g Ue -BRg™ e , for example, we minimize the Ax^ o with respect to BR^ 6 ; this 
minimization is achieved simply by replacing \x\i from the sum of IV in the definition of the 
Ax^ , and then we take 6.3 for the resulting Ay 2 (BR 6 *, BR fi *) to have the projected contour 
of90%CL. 

In Fig. [71 the region where all of the CP-conserving cases I, II, III, and IV are rejected is 
given by calculating the Ax C py function. The Ax^py is 

{/jytruc _ jyfilA 2 ] 
V i-f } . (26) 



=ee,efi,fifi 



We use Ay^py in the same way as A% 



2 

mo " 
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A. 2 Allowed region of ipi, (f2, and uiq 

In Fig. [8], the allowed region of (ip%, <p 2 , m o) for a given set of BR^ r , ue is obtained by 

AxkwW^o) = mini £ (^_zM_| . (27) 

Kll'=ee,e(j,, fj.fi " y 

iVj^ is a function of (ipx, <fi2, m o, x fit ). In order to project the 3-dimensional allowed region onto 
2-dimensional space, we minimize Ax^ llow with respect to a parameter. Then, we take 6.3 for 
the resulting A% 2 to obtain the 90%CL contour. 
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